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O I Abstract. Microrheology is a branch of rheology having the same principles as 

^ ' conventional bulk rheology, but working on micron length scales and nl volumes. 

Optical tweezers have been successfully used with Newtonian fluids for rheological 
O ' purposes such as determining fluid viscosity. Conversely, when optical tweezers are used 

to measure the viscoelastic properties of complex fluids the results are either limited 
to the material's high-frequency response, discarding important information related to 
the low-frequency behaviour, or they are supplemented by low-frequency measurements 
performed with different techniques, often without presenting an overlapping region 
of clear agreement between the sets of results. We present a simple experimental 
procedure to perform microrheological measurements over the widest frequency range 
ly-^ , possible with optical tweezers. A generalised Langevin equation is used to relate the 

f^ ' frequency-dependent moduli of the complex fluid to the time-dependent trajectory of 

^^ ■ a probe particle as it flips between two optical traps that alternately switch on and off. 
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1. Introduction 

Optical tweezers have become an increasingly important tool for the manipulation of 
micron-scale objects. Taking advantage of the optical gradient force [1] they have 
been used both as a force sensor |2l |3] and as an force actuator [H |5]. Though the 
single beam gradient trap has been around since Ashkin [6], modern holographic optical 
tweezers, which use spatial hght modulators to create multiple optical traps, are capable 
of moving multiple micron-sized objects in three dimensions [7]. Recent advances in 
SLM technology mean that optical tweezers are now capable of updating trap positions 
at hundreds of frames per second with low latency (~5ms) [8]. This can only serve to 
broaden the already extensive range of applications for which optical tweezers have been 
used, such as cell biology [9], investigations into DNA [5], and driving micro-machines 

m- 

Along with biological and physical experiments optical tweezers have also proved 
useful for mi crorheo logical experiments [H] . Microrheology is concerned with the linear 
viscoelastic response of materials at microscopic length scales. Such information is 
invaluable for the investigation of unknown biological processes as well as for increased 
understanding of the basic physics of fluids. The linear viscoelastic properties of a fluid 
can be represented by the frequency-dependent dynamic bulk modulus G*{uj). The 
dynamic bulk modulus encompasses information about both the elastic and viscous 
properties of the fluid. It is expressed in the form G*{u) = G'{u) + iG"{u) [12], where 
G'{uj) is the frequency-dependent elastic modulus and G"{uj) is the frequency-dependent 
viscous modulus [T3] . 

As is well known [HI [151 [E] , the viscosity of the medium in which a trapped particle 
is immersed is a key factor affecting the statistical distribution of the particle's time- 
dependent position. Hence the viscosity is easily determined from measurements of that 
distribution. Furthermore, surface effects such as Faxen's correction to the viscosity of a 
fluid close to an infinite surface ( [IS [161 ITTj) have also been widely studied in Newtonian 
fluids. However when microrheological measurements of viscoelastic substances have 
been made, results are either limited the high end of the frequency response [TH [191 [20] 
(omitting the low frequencies entirely) or account is taken of the low frequency response 
only by the use of complementary techniques such as rotational rheometry [21] or passive 
video particle tracking microrheology [22]. In both cases this has been done without 
showing an overlapping region, thus leaving a significant information gap. 

The aim of this paper is to present a self-consistent procedure for measuring the 
linear viscoelastic properties of materials across the widest frequency range achievable 
with optical tweezers. In particular, the procedure consists of two steps: (I) measuring 
the thermal fluctuations of a trapped bead for a sufficiently long time; (II) measuring 
the transient displacement of a bead flipping between two optical traps (spaced at flxed 
distance Dq) that alternately switch on/off at sufficiently low frequency. The analysis of 
the first step (I) provides: (a) the traps stiffness (kj, i = 1, 2) — note that this has the 
added advantage of making the method self-calibrated — and (b) the high frequency 
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viscoelastic properties of the material, to high accuracy. The second step (II) has the 
potential to provide information about the material's viscoelastic properties over a very 
wide frequency range, which is only limited (at the top end) by the acquisition rate of 
the bead position (~kH) and (at the bottom end) by the duration of the experiment. 
However, because of the finite time required by the equipment to switch on/off (i.e. tens 
of milliseconds), the material's high-frequency response can not be fully determined 
by this step. The full viscoelastic spectrum is thus resolved by combining the results 
obtained from steps (I) and (II). 

2. Analytical model 

To understand the basis of the procedure, we consider the time-dependent position f{t) 
of a bead trapped by a stationary harmonic potential of force-constant Kj. Throughout 
the first step (I) of the procedure, the bead is always found close to the centre of the 
single trap which is switched on; it makes only small deviations r(t) (of magnitude 
set by the thermal energy) away from the centre of the trap. During step (II) of the 
procedure, the detours are considerably larger, of a magnitude set by the separation Dq 
of the traps. Let us define tu to be the time at which step (II) commences, i.e. the time 
at which the traps are first switched. Subsequently, each trap remains on for a duration 
P before it is switched off and the other trap switched on. Hence the total period of 
the repeated sequence is 2P. At the instant immediately after the traps are switched 
(i.e. at t = til + nP where n = 0, 1 . . . N), the bead is typically positioned at a distance 
\f{t)\ ~ Do from the centre of the currently active trap (i.e. close to the centre of the 
trap that has just switched off). Note that coordinates are re-defined so that the bead's 
displacement r(t) is always measured with respect to the centre of whichever trap is 
currently switched on. The bead's position in three dimensions can be modelled by a 
generalized Langevin equation, 

md{t) = fnit) - [ Cit- r)v{r)dT - Kif{t), (1) 

J to 

where m is the mass of the particle, a{t) is its acceleration, v{t) its velocity and /ij(t) 
is the usual Gaussian white noise term, modelling stochastic thermal forces acting on 
the particle. The integral, which incorporates a generalized time-dependent memory 
function C(^)! represents viscoelastic drag from the fluid. 

Note that the present method does not require the two traps to be equal. They can 
be independently, but not simultaneously, calibrated during step I by appealing to the 
Principle of Equipartition of Energy, 

lksT=^-K,{r'), (2) 

where ks is Boltzmann's constant, T is absolute temperature and (r^) is the time- 
independent variance of the particle's displacement from the trap's centre. Despite the 
variety of established methods for determining the stiffness of an optical trap (e.g. using 
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the power spectrum or the drag force [TBI E])) the Equipartition method is the only 
one independent of the viscoelastic properties of the material under investigation and is 
thus essential to the calibration of a rheological measurement. 

We now examine how Eq. ([1]) evolves in the two steps (I) and (II). During step 
(I), the thermal fluctuations of the trapped bead are investigated to determine the high- 
frequency viscoelastic properties of the material through analysis of the time dependence 
of the normalised position autocorrelation function A[t) ("NPAF"): 

(r(to)f(to + r)),„ 



A{t) 



(3) 



which, by time-translation invariance, is a function only of the time interval r. Here, 
(r(to)^)ifi is the time-independent variance and the brackets {■■■)tg denote an average 
over all initial times to- 

Multiplying both the sides of Eq. ([1]) by r{to) and averaging over t^, one obtains 

-1 



A{s) 



s + 



rhj 



(4) 



ms + <^{s 

where A{s) is the Laplace transform of A{t), s is the Laplace frequency and it has been 
assumed that both the quantities {r{to)v{to))^ and (r(to)/ij(^o + t) ) vanish for all r. 



to 



We note that, for ms ^ ({s) (which is a good approximation up to MHz frequencies 
for micron sized beads with density of order Ig/cm? in aqueous suspension) and for a 
Newtonian fluid (i.e. a liquid with time- independent viscosity r/, for which Q = Gnrja), 
Eq. (jl]) recovers the well-known result for a massless particle harmonically trapped in a 
Newtonian fluid, 

A{t) ^ exp (-r,r) , (5) 

where Fj = Ki/Girari is the characteristic relaxation rate of the system and can be used 
to determine rj once both Kj and a are known. In the general case of non-Newtonian 
fluids (i.e. materials with time-dependent viscosity f]{t)) we follow Mason and Weitz 
[23] in assuming that the bulk Laplace-frequency-dependent viscosity of the fluid fi{s) 
is proportional to the microscopic memory function ({s) = 6TTafj{s), so Eq. (jl]) can be 
written as 



fjis) 



Ki 



Gna 



Ms) 



sAis] 



ms 



(6) 



Moreover, given that G*{uj) = srj{s)\^^^, the complex viscoelastic modulus G*{uj) can 
be expressed directly in terms of the time-dependent NPAF, 



G*{uj) 



Gvra 



iuAioj) 



1 — iujA{uj) 



+ 



muj 



Kj 



(7) 



where A{uj) is the Fourier transform of A{t) and, as mentioned before, the inertia term 
(mo;^) can be neglected for frequencies u ^MHz. 
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It is interesting to highlight that Eq. ([7]) is actually equivalent to Eq. (6) in 
Ref. [13]. Indeed, A{t) is directly related to the normalised mean-square displacement 

(Ar2(r)) /2 (r^) introduced in Ref. [I3| as follow: 

(Ar^(r)) {r\T)) + {r^tp)) - 2 (r-(to)f(r)) 

2 (r2) 2 (r2) 

= 1-A(r). (8) 

By performing the Fourier transform of Eq. ([8]) one obtains the relation 

Note that the quantity A{t) has the added advantage of having a well-controlled Fourier 
transform unlike the MSD{t) as discussed below. 

The second step (II) of the procedure consists of analysing the bead's transient 
displacements as it moves between two traps with separation Dq that swap their on/off 
state at times t = to + nP. Note that the duration P must exceed all of the material's 
characteristic relaxation times. We define the normalized mean position of the particle as 
D{t) = \{f{t))\ /Dq, where the brackets (...) denote the average over several independent 
measurements, but not over absolute time, since time-translation invariance is broken 
by the periodic switching. In this case. Equation ([1]) yields, in the Laplace form, an 
identical expression to Eq. (|4]) with A{s) replaced by D{s), the Laplace transform of 
D{t). Thus, as before, the complex modulus can be expressed directly in terms of the 
bead position: 



G*(uj) = -^ Y 



VlT^a 

j=l,2 



KiicoD(uj) n 

+ mu 



1 — iuIDiu 



(10) 



where D{<jj) is the Fourier transform of D{t) and the sum takes account of the linearity 
of the measurements performed with both traps. 

Note that both functions A{t) and D{t) are expected to have the limits A(0) = 
D{Qi) = 1 and ^(oo) = D{oo) = 0. This turns out to be very useful when applying the 
Fourier transform during data processing, as shown below. 

In principle, Eqs. ([7]) and ( ITOj) are two simple expressions relating the material's 
complex modulus G*{uj) to the observed time-dependent bead trajectory r(t) via the 
Fourier transform of either f{t) itself (in Eq. ( ITOl) ) or the related NPAF (in Eq. ((Tj)). In 
practice, the evaluation of these Fourier transforms, given only a finite set of data points 
over a finite time domain, is non-trivial since interpolation and extrapolation from those 
data can yield serious artefacts if handled carelessly. 

The N experimental data points [tk, A{tk)) or (t^, D{tk)), where A; = 1 . . . A^, extend 
over a finite range of times t, exist only for positive t and need not be equally spaced. In 
order to express the Fourier transforms in Eqs. ([7]) and fITOl) in terms of these data, we 
adopt the analytical method introduced in Ref. [21]. In particular, we refer to Eq. (10) 
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of Ref. [23] which is equally applicable to find the Fourier transform g{u)) of any time- 
dependent quantity g{t) sampled at a finite set of data points {tk,gk), giving 

(^i-^(O)) 



u'g{u)=zug{0)+{l-e-'-'^) 



ti 



N 



where 5^00 is the gradient of g(t) extrapolated to infinite time. Also g{0) is the value 
of g{t) extrapolated to t = 0^. Identical formulas can be written for both A{u) and 
D{uj), with g replaced by A and D respectively. It is a strength of our new procedure 
that both of the extrapolated quantities cjoo and g{0~^) (which assumed the role of fitting 
parameters in Ref. ^4]) in this case assume known values (5^00 = and g{0~^) = 1) given 
by the limits of the functions A{t) and D{t). Like the method presented in Ref. |24j . 
the present procedure also has the advantage of removing the need for Laplace/inverse- 
Laplace transformations of experimental data [25] . 

3. Experimental Setup 

We have validated the experimental procedure as described by Eqs. (jTj) and (fTOl) . via 
Eq. (ITTi) . by measuring both the viscosity of water and the viscoelastic properties of 
water-based solutions of polyacrylamide (PAM, fiexible polyelectrolytes, M^ = 5 — 
6 X 10^ g/mol, Polysciences Inc.) using optical tweezers as described below and shown 
schematically in Fig. [H 

Trapping is achieved using a CW Ti:sapphire laser system (M Squared, SolsTiS) 
which provides up to 1 W at 830 nm. Holographic optical traps are created via the use 
of a spatial light modulator (Boulder XY series) |8J in the Fourier plane of the optical 
traps. The tweezers are based around an inverted microscope, where the same objective 
lens, 100 X 1.3NA, (Zeiss, Plan-Neofiuor) is used both to focus the trapping beam and 
to image the resulting motion of the particles. Samples are mounted in a motorized 
microscope stage (ASI, MS-2000). Particles are imaged using bright-field illumination. 
We use a Prosilica GC640M camera to view the trapped particle and our own suite 
of camera analysis software written in Lab VIEW (available from j26]) to measure the 
position of the trapped particle in real time at up to 2 kHz frame rate [27] . 

4. Results 

The normalized position autocorrelation functions A{t), measured from the stochastic 
fiuctuations of beads optically trapped in two different fluids are shown in Fig. El 
together with the function predicted for a simple Newtonian fluid. In the Newtonian 
case, it is expected that A{t) decays as a single exponential with a characteristic 
relaxation rate related to the trap strength, bead size and fluid viscosity i.e. Eq. (|5]). 
The agreement between the data and prediction for water is good. On the other hand. 
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Figure 1. Experimental setup. From left to right. (LZ) Tiisapphire laser system, (LI 
and L2) Beam Telescope, (Ml and M2) Folding Mirrors, (SLM)Boulder Fast SLM,(L3), 
(PBS) Polarising beam splitter cube, (M3) Mirror, (CAM) Prosilica Fast camera, (O) 
Objective lens, (CL) Condensing optics, (B) 250W halogen bulb. 



in the case of a non- Newtonian fluid, where the viscosity is time- dependent, it is not 
guaranteed that Eq. @ could be resolved (i.e. inverse-Laplace transformed) into a simple 
form like Eq. ([5]). However, this is no hindrance since the viscoelastic moduli will be 
found via the analysis of the normalized position autocorrelation function. In Fig. |3]we 
compare the impulse response (i.e. step II of the procedure) of a 5/im diameter bead 
suspended in water and in an aqueous solution of PAM at 1% w/w (a non- Newtonian 
fluid), with a duration between flips of P = 20 s and a trap centre-to-centre separation 
of Dq = 1.6fim, giving Do/a = 0.64. In order to guarantee the linearity of the confining 
forces exerted by the two optical traps, the distance between them was always chosen to 
be no more than 80% of the bead radius, Dq < 0.8a p. Although Brownian statistical 
fluctuations appear in the figure, the difference between the viscoelastic natures of the 
two fluids is clear. Indeed, while a bead suspended in water flips from one trap to the 
other almost instantaneously, the same bead in the PAM solution takes much longer 
(a few seconds) to flip. In order both to evaluate D(t) = \{r(t))\ /Dq and to reduce 
noise caused by Brownian fluctuations, the transient measurements were averaged over 
twenty flips, with the resulting curves shown in Figure HI Experimentally, the switching 
process of the two traps is controlled by means of a spatial light modulator (SLM) which 
alternately creates an optical trap in one of two positions. It is important at this point 
to note that there is a short but finite time for which both traps exist simultaneously [8j, 
due to the finite time required by the SLM's display to update the holographic pattern. 
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This makes the switching process not exactly binary. However, this will of course only 
affect the high-frequency results obtained during this step which will ultimately be 
neglected in favour of the high-frequency response from step (I). 

Wideband microrheological measurement are obtained from the optical tweezers 
by combining the frequency responses obtained from both steps (I) and (II) of the 
procedure. In particular, the material's high-frequency response is determined by 
applying Eq. ([7]) {via Eq. f fTTj) with A/^ replacing g^) to the A{t) measurements (in 
which low-frequency information tends to be very noisy); whereas, the low-frequency 
response is resolved by applying Eq. ( !T0|) {via Eq. (TTT]) with D^. replacing gk) to the 
data describing the bead's transient response to the flipping traps (in which the high- 
frequency response is limited by the performance of the SLM). 

Typical results for both Newtonian and non-Newtonian fluids are shown in Figures [5] 
and El respectively. In both the cases, it is evident that, although there is some noise 
in the frequency domain which comes from genuine experimental noise in the time- 
domain data, there is a clear overlapping region of agreement between the two methods 
which makes the whole procedure self-consistent. Moreover, it confirms the ease with 
which the low-frequency material response can be explored right down to the terminal 
region (where G' oc u"^ and G" oc u). This is the current limitation for microrheological 
measurements. 

5. Conclusions 

In summary, we have presented a self-consistent and simple experimental procedure, 
coupled with a data analysis methods, for determining the wide-band viscoelastic 
properties of complex fluids using optical tweezers. This method extends the range of the 
frequencies previously available to optical tweezers measurements. In fact, the accessible 
frequency range is limited only by the experiment length and by the maximum data 
acquisition speed (10s of MHz for a quadrant photo-diode). This allows access to the 
material's terminal region enabling micro-rheological measurements to be performed on 
complex fluids with very long relaxation times, such as those exhibiting soft glassy 
rheology |28]. The method provides a simple yet concrete basis on which future 
viscoelastic measurements may be made on both biological and non-biological systems 
using optical tweezers. 
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Figure 2. The normalised position autocorrelation function vs. lag-time of a 5/iTO 
diameter bead (squares) in water (with k — 2.7^N/m) and (circles) in a water-based 
solution of PAM at concentrations of 1 % w/w (with k = 2.2^N/m). The continuous 
and dotted lines represent Eq. ([5]) for a bfim diameter bead in water at T = 25°C with 
n — 2.7iiN/m and k — 2.2iiN/m, respectively. 
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Figure 3. The trajectory of a 5/zm diameter bead flipping between two optical traps 
Ki (bottom) and K2 (top) repeatedly switching after a duration P = 20 s. The bead 
is suspended in (squares) water (with ki = 2.7 and K2 — 2.biiN/m) and (circles) 
a water-based solution of PAM at concentrations of 1 % w/w (with ki = 2.1 and 
K.2 = 2.2fj,N/m). 
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Figure 4. The normilized mean position of all step-down data shown in Fig. [31 i.e. 
when simultaneously trap 2 (top) switches off and trap 1 (bottom) switches on. 
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Figure 5. Storage (G") and loss (G") moduh of water vs. frequency, analysed using 
both Eq. ([7|) (high frequencies) and Eq. ([TU]) (low frequencies) applied directly to the 
experimental data presented in Fig. [2] and Fig. |4l respectively. The lines represents 
the expected limiting behaviour of the moduli when the material reaches the terminal 
region: G' oc w^ and G" oc u. 
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Figure 6. Storage (G") and loss (G") moduli us. frequency of a solution of 1% w/w of 
PAM in water measured by means of both Eq. Q (high frequencies) and Eq. (fTO|) (low 
frequencies) applied directly to the experimental data presented in Fig. [2] and Fig. SI 
respectively. The lines represents the expected limiting behaviour of the moduli when 
the material reaches the terminal region: G oc w^ and G" oc lo. 
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